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How did interpersonal violence develop in early human societies? Given 

that homicide records are only available for the more recent period, much of 
human history remains outside our purview. In this paper, we study violence 
trends in the very long run by exploiting anew dataset on cranial trauma and 


weapon-related wounds from skeletons excavated across the Middle East, 
spanning the pre-Classical period (around 12,000-400 BcE). The dataset 
includes more than 3,500 individuals. We find evidence that interpersonal 
violence peaked during the Chalcolithic period (around 4,500-3,300 BCE). 
It then steadily declined during the Early and Middle Bronze Ages (around 
3,300-1,500 BCE) and increased again between the Late Bronze and the 

Iron Age (1500-400 BCE). By documenting variations in violence patterns 
across a vast temporal and geographical scale in an incredibly rich historical 
setting, we broaden perspectives on the early history of human conflict. 


How did violence develop in the very long run? Some recent studies 
argued for asecular trend toward less violence over centuries or even 
millennia’ °. However, the prevalent claim that interpersonal violence 
(assault, murder, slavery, torture, despotism, cruel punishment, violent 
feuds and so on) has declined over millennia—and more rapidly since 
the Enlightenment—has been variously challenged*». But so far, no con- 
sensus has been reached. In fact, because homicide statistics are only 
available for recent periods** and because conflict records are subject 
to reporting biases‘, our understanding of long-term violence trends 
is still limited, the further back we go in time. This situation has gener- 
ated highly polarized narratives about the early history of violence”. 
Indeed, the scholarship on the origins of warfare split into two major 
factions: the so-called ‘pacified past’ group or ‘doves’—thinking that 
violence in human history picked up only after the Agricultural Revolu- 
tion and the rise of states; and the other viewing warfare as pervasive 
in prehistory, that is the ‘hawks? ™. Although the hawks/doves division 
fuelled much debate between the 1990s and early 2000s, nowadays a 
more nuanced picture is emerging. 

Lately, social scientists, anthropologists, biologists and archaeolo- 
gists interested in historical conflict have started to converge in terms 
of methods and scopes. Studies investigating the causal mechanics of 
violence, often adopting a cross-cultural approach, have started to bet- 
ter document variations in conflict patterns through the accumulation 


of new skeletal information’* ®. Recent progress in forensic investi- 
gation techniques applied to human skeletal remains prompted the 
proliferation of standardized databases of excavated skeletal data? ”, 
contributing tremendously to studying variations in human violence 
across different spatial, social and temporal scales, which is crucial 
to correctly frame the role of conflict in human evolution”. 

We build upon these theoretical and methodological advances to 
construct and analyse a dataset of skeletal evidence which allows us to 
study long-term violence trends in the ancient Middle East, spanning 
almost twelve millennia and covering all its major regions (Fig. 1). Our 
evidence-based contribution champions the potential for biocultural 
interpretation of large bioarchaeological datasets when answering 
bigger questions about past violence. Understanding differences in 
levels of violence in the past is not only a central issue in social 
sciences!” but also in anthropology, archaeology and history”?”?”*. 

We formulate several proposals, on the basis of the existing lite- 
rature as follows. (1) During the Neolithic period, the population den- 
sity was still very low and therefore potentially conflicting groups 
could easily avoid each other. During this period, we would expect 
relatively low incidence of violence, roughly in the range or slightly 
below the European violence levels observed by ref. 25. (2) During the 
Chalcolithic period, when the first cities emerged, we would expect 
violence to be more common. This was a period that witnessed stressful 
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Fig. 1| Distribution of Middle Eastern archaeological sites that provided bioarchaeological evidence for the present study. In the map, all individual sites are 
shown, plus the four regions of Turkey, Mesopotamia, Levant and Iran. Base map adapted from ref. 98, CC BY 4.0 (https://www.gmrt.org/about/terms_of_use.php). 


changes such as urbanization and the emergence of protostates”’. 


(3) During the Early and Middle Bronze Age, the synergetic effect of 
state, law and trade development reduced interpersonal violence. 
(4) During the Late Bronze and Iron Ages, we would expect higher vio- 
lence due to the frequently postulated ‘collapse’ phenomena and the 
only gradual recovery of conflict-reducing state activities. 


Results 

Scope of the study 

Here, we report the findings of a study based on the analysis of anew 
dataset on violence-related injuries in seven Middle Eastern countries 
(Turkey, Iraq, Iran, Syria, Lebanon, Israel and Jordan). We study no less 
than 3,539 well-provenanced skeletons dated to the period from the 
Mesolithic to the onset of Classical Antiquity, around 12,000-400 BCE 
(Table 1)”””’. We define the intensity of interpersonal violence as the 
share of individuals exhibiting weapon-related wounds or cranial trau- 
mas (Supplementary Methods). We combine primary data analyses 
andthe re-analysis of published material that used different reporting 
systems (ways of representing the observations using numerical val- 
ues), as well as published sources of cranial trauma and weapon wound 
aggregates that do not always provide a catalogue of raw data. We need 
to emphasize that perfect standardization is impossible, especially with 
the last type of source (for a detailed discussion of the limitations, for 
example arising from the crude-versus-‘true’ rates, see Supplementary 
Methods). However, we carefully assess the robustness and plausibility 
of our findings and observe them to be informative. 

The pre-Classical Middle East is a key laboratory for our exercise 
because it allows us to study changes in violence trends across several 
crucial historical transformations, such as the transition from forag- 
ing to farming, the Urban Revolution, the rise of agrarian states and 
their evolution up to universalistic empires, globalism and so on*°; 
and also because no systematic studies of comparable scale have been 
attempted before in this setting. 

Our results inform both the social scientific literature seeing the 
emergence of early states as a social technology instrumental for 
reducing violence” and the debate about the origins of violence in 
archaeology and anthropology, where an influential view is that increas- 
ingly complex political organization brought about higher intensity of 
warfare and coercive capabilities, with higher casualties and burdens 
for the population”. Our findings provide three advances over the 


extant literature. First, our approach combines methods in historical 
econometrics and anthropology for reliably assessing violence lev- 
els across the full spectrum of the historical experience—provided 
that sufficient bioarchaeological data are available—thus enhancing 
replicability. Second, we fill a gap in scholarship by shedding light on 
12,000 years of conflict in an exceptionally important historical setting, 
the pre-Classical Middle East. Third, we document that, although the 
rise of first states was accompanied by a dramatic increase in strife, 
the development of ‘mature’ Bronze Age states went in parallel with a 
strong reduction of interpersonal violence, which aligns with the pro- 
posal that by building state capacity, early institutions could provide 
some key benefits to their populations. 


Ashort history of the Ancient Middle East 

We begin with a short historical overview. How did some core aspects 
of ancient Middle Eastern societies develop over the time frame from 
12,000 to 400 BCE? We define these ‘core aspects’ as (1) the ability of 
the state to constrain interpersonal violence, (2) population density, 
(3) climate (with its impact on harvests) and (4) military development. 

During the Neolithic period, population density remained low, 
despite the gradual diffusion of farming technologies”. However, 
during the Late Chalcolithic period (from around 4,400 BCE), a trend 
towards drier conditions was one of the factors potentially triggering 
urbanization, with population flowing from the countryside to major 
sites rapidly turning into cities, where emerging elites started to com- 
pete over the control of resources’. Crowding and rising inequality 
may have triggered conflicts following appropriation of other groups’ 
resources resulting in higher interpersonal violence”. The Chalcolithic 
period saw also the rise of the first protostates*. Following the proposal 
that these were low-capacity states””*, their ability to settle conflicts 
within their populations was seemingly still quite limited. 

The following Early Bronze Age and Middle Bronze Age (hereafter, 
EBAand MBA; around3,100-1,550 BCE) were the first ‘Golden Age’ periods 
in Mesopotamia, Levant and partially in the neighbouring regions. 
City-states and empires developed increasing state capacity which we 
define as the ability of the state to reach its aims via taxation, constrain- 
ing interpersonal violence and other classical roles of the state. This 
is attested by increasingly ambitious building programmes (fortifica- 
tions, temples and palaces*”*°) and steady territorial expansion*’. 
Furthermore, arelatively high sophisticationin settling conflicts using 
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Table 1| Number of cases for violence shares by region and 
period 


Region Neolithic Chalcolithic EBA MBA LBA IA Total 
Iran 15 127 (2) 209 353 
Levant 37 75 29733 (5) 73 520 
Mesopotamia 56 (12) 90 66 51 51 326 
Turkey 1,078 367 786 100 (9) 2,340 
Middle East 1186 581 1175 199 65 333 3,539 


Regions: Iran; Levant (West Syria, Lebanon, Israel and Jordan); Mesopotamia (Iraq and East 
Syria); Turkey; Middle East, sums. Periods (all BCE): Mesolithic/Neolithic, 12,.000-4,500; 
Chalcolithic, 4,500-3,300; EBA, 3,300-2,000; MBA, 2,000-1,550; LBA, 1,550-1,200; IA, 
1,200-590. Brackets indicate low number of underlying cases. 


ideological, bureaucratic and legal means was achieved. Codified legal 
systems”, ramified bureaucracies and strongly centralized military 
control“, the diffusion of intercity religious festivals***°, cooperative 
projects” and risk-sharing activities contributed to transforming fol- 
lowers into citizens**”’, which was instrumental for controlling social 
violence. This is further attested by the fact that, although the transi- 
tion between the EBA and MBA was characterized by high population 
thresholds”, extended drought” (the BP event of 4.2 thousand years 
ago (ka)) and by a take-off of military activity (in terms of intensity 
and development of bronze weaponry”), interpersonal violence 
kept decreasing. 

Moreover, during the EBA-MBA, long-distance trade flourished**, 
as merchants were relatively protected against expropriation by 
rulers—especially private merchants in Mesopotamia”. This prosperity 
could be sustained for several centuries, even in the face of increasingly 
arid conditions” and rising inequality°®®. In fact, the combination 
of declining crop yields, increasing privatization and monetization 
of economic activities, impoverished farmers, while private traders 
and businessmen were able to amass increasingly large economic and 
political power due to their roles as money-lenders, tax-farmers and 
enterpreneurs”. As to the military sector, the eclipse of the regional 
Akkadian and Ur III empires was followed by a balkanized landscape of 
city and territorial states that competed for supremacy, often through 
warfare“, Noteworthy is the professionalization of army service 
in this phase and the increasing reliance on mercenary troops for 
fighting wars”, which may have contributed to lower the occurrence 
of violence-related injuries in the social sectors not directly involved 
in warfare”. 

Finally, the transition from the Late Bronze Age to the Early Iron 
Age (hereafter LBA and IA; around 1,550-1,100 BCE) saw a temporary 
flourishing of the Levant region with the rise of the Hittite and Mitannian 
powers”, followed by a substantial crisis of the Bronze Age societies. 
This was propelled by a 300-years drought, climate-induced migra- 
tory movements and general economic decay®* ©. Many famous cities 
were burnt by foreign invaders® and climatic worsening reduced the 
population density strongly, especially in Mesopotamia®. The spread 
of iron weapon technology and mounted warfare revolutionized war- 
fare through the creation of iron-armed and armoured shock troops, 
leading to a change in social structure; hegemonic empires expanded 
inan unprecedented way” and this also triggered a proper ‘arms race’ 
in Eurasian societies’. Consequently, in the ‘Age of Empires’, a larger 
share of the population participated—directly as mass-mobilized sol- 
diers or indirectly as civilian casualties or deportees—in warfare”. 
Another byproduct of the formation of empires controlling vast terri- 
tories and endowed with unparalleled powers of extraction and mobi- 
lization was a greater and even more unequal wealth accumulation, 
especially concentrated in very large imperial centres’. Repressive 
policies, large-scale warfare and unprecedented levels of inequality 
may also have had implications for interpersonal violence. Widespread 
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Fig. 2| Trends in violent trauma in the study sample. Regions shown are: Iran 
(n= 353); Levant (n = 520, West Syria, Lebanon, Israel and Jordan); Mesopotamia 
(n= 326, Iraq and East Syria); Turkey (n = 2,340); the Middle East, average. Periods 
(all BCE) are: Neolithic (including Mesolithic), 12,.000-4,500; Chalcolithic, 
4,500-3,300; Early Bronze, 3,300-2,000; Middle Bronze, 2,000-1,550; Late 
Bronze, 1,550-1,200; Iron, 1,200-400. We justify the combination of the two 
periods of the Mesolithic and the Neolithic by the fact that the violence rate 

was almost identical for the two. Dropping the Mesolithic would not change 

the results substantially. For Late Bronze Turkey, the underlying number of 
cases is 9, all other cases represent more than 30 observations (the other low- 
number observations in brackets in Table 1 were excluded). We also performed 
arobustness test of this figure in which we include or exclude cases with low 
numbers of observations. We discuss the results in the Supplementary section 
on robustness and sensitivity tests. We performed a number of two-sided t-tests 
toassess the statistical significance, for example, between the violence in the 
Neolithic and the Chalcolithic, using the number of underlying individuals 

as weights. For the differences between Neolithic and Chalcolithic for 
Mesopotamia, Iran, Turkey and the whole Middle East: P < 0.001. Chalcolithic 
versus Early Bronze Age: for Levant P< 0.001, for Iran and Turkey (small numbers 
of observations)—Iran P = 0.052, Turkey P= 0.637; for the Middle East P< 0.001. 
Early versus Middle Bronze Age: Levant P < 0.001, Turkey P= 0.789; the whole 
Middle East P = 0.224. Middle versus Late Bronze Age: Turkey P< 0.001; Middle 
East P= 0.042; the others could not be calculated for these periods. Late Bronze 
versus Iron Age: Mesopotamia P < 0.001; the whole Middle East P = 0.482 
(Supplementary Table 1). 


resistance to imperial power in the form of rebellions and insurgencies 
is attested inthe written sources of the time, especially in the peripheral 
and remote areas of the empire, as well as intra-elite fighting”. 


Empirical findings 

We assigned each observation to one of the six time periods, namely 
the Mesolithic/Neolithic period (from around 12,000-4,500 BCE), the 
Chalcolithic (around 4,500-3,300), the EBA (around 3,300-2,000), the 
MBA (2,000-1,550), the LBA (1,550-1,200) and finally the IA (1,200- 
400). The Mesolithic/Neolithic period is dominated by the much larger 
Neolithic sample; hence we summarize it under ‘Neolithic’ in the fol- 
lowing (Fig. 2). We follow the archaeological periodization inline with 
previous influential studies*’. We define ‘Levant’ as the rain-fed region 
encompassing modern western and northern Syria, Israel and Lebanon; 
whereas ‘Mesopotamia’ corresponds to the region where irrigated 
farming was practiced, that is modern Iraq and eastern Syria, up to 
the area of Tell Ashara (Euphrates/Habur confluence; Fig. 1). Ancient 
Turkey corresponds to modern-day Turkey, as does ancient Iran to the 
modern country. 

The empirical results suggest that there was a high level of 
violence in the Chalcolithic period and elevated levels during the Iron 
Age (Fig. 2). Instatistical tests, we observe that the differences between 
the Neolithic and Chalcolithic are statistically significant, with a Cohen’s 
d effect size of 1.193 (b is the coefficient, z-statistic with degrees of 
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freedom in brackets, Cohen’s d effect size calculated as coefficient 
relative to standard deviation, confidence intervals (CI); we consider 
model 3 as normality can be confirmed; see notes to Table 2). Neolithic 
versus Chalcolithic: b = -1.332, z(70) =—-4.33, P= 0.000, Cohen’s d effect 
size 1.193, 95% CI -1.934, -0.729 (Table 2, column 3). Similarly, the dif- 
ferences between Chalcolithic and EBA-MBA were statistically signifi- 
cant. Notably, the data from the EBA-MBA show a consistently lower 
level of violence in the region than before (Early Bronze b = —0.800, 
2(70) = -3.22, P= 0.001, Cohen’s d effect size 0.741, 95% Cl -1.287, -0.312; 
Mid-Bronze b =-1.815, z(70) = -2.55, P= 0.011, Cohen’s d effect size 
2.852, 95% CI -3.210, -0.420). Moreover, the Iron Age was different from 
the Chalcolithic: b =-1.114, z(70) = -4.38, P= 0.000, Cohen’s d effect 
size 0.997, 95% Cl -1.612, -0.615. We observe significant difference for 
burial contexts that were related to battles (b = 1.157, z(70) = 2.600, 
P=0.009, Cohen’s d effect size 1.134, 95% CI 0.283, 2.031) and royal 
context (b = 1.861, z(70) = 3.970, P= 0.000, Cohen’s d effect size 1.824, 
95% CI 0.941, 2.781). The coefficients are reported in logs (on the loga- 
rithmic transformation; Methods). If we want to express the effect as per 
cent, we need to calculate the exponential. For example, the coefficient 
of MBA (relative to Chalcolithic) is exponential(1.815) = 6.141, witha 
negative sign, more than 6% decline in violence. 

As to regional distribution, we detect an occurrence of traumas of 
8.4% in Turkey, 18.6% in Iran and 10.8% in Levant; and for Mesopotamia, 
6.3% (Fig. 3). However, the error bars (representing standard errors) 
indicate that the data for each region vary quite strongly. One poten- 
tial reason is that changes over time were quite substantial (as we see 
below in the data) and this variation enters the total standard error of 
each region. 


Discussion 

In the following, we interpret these changes on the basis of our ini- 
tially formulated proposals. The elevated levels during the Chalcolithic 
and IA periods are not unexpected. Indeed, the Chalcolithic period 
witnessed the shift from occasional feuding to large-scale organized 
conflict which has been connected to the emergence of the first centra- 
lized protostates”*”°”. Similarly, the transition to the IA experienced 
major upheavals—a 300-year drought, population displacements, 
resource stress and state dissolutions, especially in the Levantine 
regions and Iran°*—followed by the rise of superpowers, such as the 
Assyrian empire, that expanded through military campaigns, forced 
taxation and mass-scale deportations~”’. A possibly similar mechanism 
is identified by ref. 15 for the Inca period Nasca highlands in Peru, where 
a reduction in resource availability due to negative climate shocks 
triggered competition and escalation in lethal conflict in populations 
with high demographic thresholds. 

However, a substantial and continuous decline in lethal aggres- 
sions was experienced between 3,000 and 1,500 BCE, during the EBA 
and especially MBA. Notably, despite state formation/decline episodes, 
droughts and warfare**—especially across the transition between the 
third and second millennia, when a shift towards drier conditions 
increased pressure on resources and was one of the triggers of demo- 
graphic changes”’’”—the Bronze Age Mesopotamian societies managed 
to reduce internal violence. They partly redirected conflict towards 
external enemies, as indicated by lower levels of trauma in Mesopota- 
mia and by the match between the concentration of traumas and the 
royal inscriptions documenting the warfare campaigns of Mesopota- 
mian kings along the eastern borders”*”. 

An influential literature on historical violence trajectories 
posits that the combined effect of the development of state capacity 
and trading systems is possibly the main force behind the declining 
homicide rates in Medieval Europe. Three key historical facts from 
our sample seem to support this interpretation (see also Supple- 
mentary Appendix A). First, although protostates emerged very 
early”, itis only around 2,500 BCE that Mesopotamian polities morph 
into territorially centralized kingdoms. Later, they developed into 


1,3,78 


empires with extensive administrative apparatuses and strong polic- 
ing capabilities*’”’. Second, written law codes and legal remedies 
for private rights—enforced by formal courts—expanded the most 
between 2,400 and 1,700 BCE (ref. 80). Third, empirical analysis of 
long-distance exchanges indicates that trade expansion peaked 
between the late EBA and the early MBA**>**". Notably, together with 
the post-2,000 BCE ‘trade revolution’, cuneiform writing expanded 
dramatically both geographically and socially, indicating rising levels 
of literacy®. These changes were flanked by the increasing profession- 
alization of Mesopotamian armies”, which possibly contributed to 
confining interpersonal violence to specific social arenas. Thus, state 
consolidation and the expansion of legal and commercial systems 
may have exerted an effect of lessening the propensity for violence 
inthe sample. 

Moreover, we assessed whether the changes towards less inter- 
personal violence during the EBA-MBA might have been rather cor- 
related with inequality, climate, population fluctuations or major 
military innovations. Consistent with our main proposal, we found that 
these alternative explanations do not seem to bear much weight. First, 
inequality—a possible driver of violence®’—steadily increased during 
the Bronze Age” and it is therefore at odds with our findings. Second, 
when we formally assess the relationship with droughts, demography 
and the technological innovations (see section on Main results of 
regression analysis below) associated with the ‘cavalry revolution” 
(Table 2), we did not observe a strong correlation with lethal aggression. 
Overall, this test, although limited to the violence covariates that are 
most efficiently measurable in archaeological and historical records, 
suggests that it is difficult to envision mechanisms different from the 
ones we take into consideration. 

Similarly, and compatible with our regional results, others find 
evidence of declining violence in the most strongly organized Pueblo 
societies (that is, when their state became more capable), inthe south- 
western United States, that could be linked to the establishment of 
intervillage confederacies, whereas the expansion of aggressive polities 
and societal collapse seemingly increased strife’*”*. 

Although we cannot draw any causal conclusion or generalization 
from the observed synchronisms, we can confirm for this specific case 
of the Ancient Middle East that the violence decline took place at atime 
whenearly states achieved substantial capacities to reduce conflictsin 
their societies*°*' and then increased again between the LBA and the 
IA (1,500-400 BCE), when major crises and the expansion of empires 
caused upheavals across the whole region”. Given that the number 
of observations is quite low for the LBA, we need to take the results for 
this period with caution (the statistical significance of the LBA is also 
not givenin columns 2 to7 of Table 2). The increase in violence between 
the MBA andthe LBA is plausible, if we see the higher violence level of 
the LBA and the IA together. The early IA especially shares character- 
istics with the LBA. 

Through a standardized methodology for inferring cranial 
and weapon-related trauma in sufficiently well-preserved skeletons 
(Supplementary Methods), we find that interpersonal violence peaked 
during the Chalcolithic period (around 4,500-3,300 BCE). Then violence 
steadily declined during the EBA-MBA (around 3,300-1,500 BCE). 
In an additional analysis (Supplementary section on robustness 
and sensitivity tests), we observe that the Chalcolithic increase or the 
EBA-MBA decline of violence was not only a proto-urban phenomenon, 
as the splitting into urban and rural effects results in overall significant 
effects but not a separate urban effect. 

We also assess (and reject) the possibility below that ‘warrior 
idealization’ might have biased the trends. This idealization might be 
reflected ina high number of bronze weapons which increased during 
the Middle Bronze Age (Fig. 4; Methods)’. Given that this was a period 
of declining personal violence, a ‘warrior idealization bias’ seems not 
to have implied a higher rate of cranial trauma and weapon wounds 
from violent conflicts. 
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Table 2 | Regression analyses of violence in the periods studied using coefficients of variables potentially related to 


violence 
Model (1) (2) (3) (4) (5) (6) (7) 
Neolithic -1.951** -1.779** -1.332** -1.315** -1.303** -1.293** -1.297** 
(0) (0) (0) (0) (0) (0) (0) 
Early Bronze Age -1.398** -1.213** -0.800** -0.752** -0.739* -0.737* -0.755** 
(0) (0) (0.001) (0.007) (0.011) (0.010) (0.007) 
Middle Bronze Age -1.730 -1.417 -1.815* -1.714* -1.935* -1.947* -1.881* 
(0.051) (0.075) (0.011) (0.018) (0.013) (0.015) (0.027) 
Late Bronze Age -0.704* 0.095 0.715 0.738 0.666 0.665 0.726 
(0.040) (0.867) (0.377) (0.366) (0.430) (0.432) (0.442) 
Iron Age -0.821* -1.050** -1.114** -0.867** -0.759* -0.733* -0.753 
(0.044) (0) (0) (0.002) (0.035) (0.042) (0.061) 
Battle sites 1157** 1.260** 1.248** 1.249** 1.182* 
(0.009) (0.006) (0.007) (0.007) (0.011) 
Royal interments 1.861** 1.610** 1.601** 1.591** 1.674** 
(0) (0.001) (0.001) (0.001) (0.005) 
Urban 0.551 0.568 0.524 0.533 0.505 
(0.199) (0.219) (0.270) (0.285) (0.300) 
Climate shock -0.728 -0.735 -0.752 -0.702 
(0.120) (0.117) (0.097) (0.166) 
Iron 0.386 0.400 0.356 
(0.463) (0.443) (0.536) 
Horse -0.062 -0.093 
(0.861) (0.781) 
Population size 0.008 
(0.604) 
Regional fixed effects No Yes Yes Yes Yes Yes Yes 
Constant -1.829** -1.297** -1.939** -1.798** -1.943** 19654 -1.996** 
(0) (0) (0) (0) (0) (0.001) (0.001) 
Observations 82 82 82 82 82 82 82 
P value (normality) 0.001 0.023 0.369 0.517 0.824 0.824 0.756 


* and ** indicate 5% and 1% significance levels, respectively, using Huber-White-Sandwich standard errors (hence heteroskedasticity is not an issue). Clustering at region/period level. All P 
values are based on two-sided tests. Seven multiple Tobit regressions show coefficients of various variables potentially related to violence. The first line shows seven different coefficients 

for the Neolithic (versus Chalcolithic), taking into account other combinations of potential correlates. Potential correlates are listed on the left. If a number and bracket term is in columns 1 to 
7, the variable is included in the specific model. Normality of the residuals is given for models 3 to 7 (not 1 and 2); we do not interpret the significance in columns 1 and 2. We use logarithmic 
transformation; following the standard procedure, we added a small constant (0.05; using other small numbers does not change the results notably). Ordinary-least-square estimation yielded 
almost the same results. Heteroskedasticity is not an issue but we tested equality of variances and found no violation of this assumption. A multiple comparison adjustment is not used, as we 
test only one dependent variable and a small set of theoretically justified covariates (for Holm-Bonferroni adjustments, see Supplementary Appendix C, showing almost the same results). 
The constant always refers to the Chalcolithic period. In column 2, the model considers region-specific effects. In column 3, also royal interments, battle site and urban/rural differences 

are accounted for (the constant refers to rural places). In columns 4 to 7, climate shocks, iron, horse and population size are added. The variables battle, royal, urban, iron and horse are 1, if 
their specific characteristic applies (being a battle site and so on) and O otherwise. Climate shock is defined as the share of the period a specific place was affected by droughts and lower 


temperatures; ‘size’ is estimates of the population size of a specific site during a specific period. 


These findings are not only useful for providing an empirical 
grounding to the history of violence in the pre-Classical Middle East 
but they also inform debates in anthropology, archaeology and social 
sciences about the nexus between conflict and premodern state forma- 
tion. In particular, the emergence of early states in the ancient world 
is usually seen as a process in which political leaders, acting as violent 
entrepreneurs, establish ‘armed peace’ via coercion and absolutism’. 
Given the lack of violence time series for much of early human history 
(except for Europe in the Neolithic”), this proposal has not been sub- 
jected to formal scrutiny. By exploiting bioarchaeological evidence 
not previously compiled to quantify levels of interpersonal violence 
in the deep past, we document more accurately when and in what 
circumstances ancient Middle Eastern societies became effective in 


controlling interpersonal violence. We find that, although the emer- 
gence of protostates coincided with a substantial violent spike in our 
sample, the timing of violent rates reduction is compatible with the 
theory that once early states achieved core capabilities—that is, terri- 
torial control, centralization of political violence and full-fledged 
legal systems—and actively supported the expansion of trade networks, 
levels of lethal violence markedly dropped (proposals 1-3). However, 
these security gains—seemingly aided by climate and epidemiological 
shocks temporarily reducing population pressure—were short-lived 
and the region witnessed major disruptions across the transition 
between the Bronze and Iron Ages, possibly because of the rise of 
more extractive powers, more widespread access to cheaper weaponry 
and mounting pressure over diminishing resources (proposal 4). 
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Fig. 3 | Share of violence in the four regions (percentage of individual human 
remains). The error bars represent standard errors of cranial trauma and weapon 
wounds, by region. The bars represent averages of violence by region. 


We conclude that the growing availability of anthropological/ 
archaeological proxies and the convergence of the application of 
statistical methods and biocultural interpretations provide new path- 
ways for interdisciplinary cooperation between social sciences and 
historical disciplines on crucial debates, such as the role of violence 
in human evolution. 


Methods 

Measuring interpersonal violence 

To study variations in violence levels in our sample, we measure the 
intensity of interpersonal violence as the share of individuals exhibiting 
weapon-related wounds or cranial traumas (Supplementary Methods; 
refs. 78,84). The use of this indicator is motivated by the fact that in his- 
torical populations, cranial traumas and weapon injuries were often the 
results of interpersonal violence, as demonstrated by the systematic 
analysis of such indicators for a skeletal sample from European regions 
over 2,000 years”. To construct our database, we collected evidence of 
antemortem cranial traumas, excluding postmortem damage. We also 
considered perimortem trauma if it could be clearly distinguished from 
postmortem damage (following ref. 78). We excluded postmortem 
damage, as postmortem fractures can be identified—again with some 
measurement error—by the amount of collagen of the bone that was 
already lost at the time of the fracture. These fractures, as far as pos- 
sible, tend to propagate at right angles to the bone surface, rather than 
an oblique angle as in the case of perimortem trauma. Often, a colour 
difference relative to the undamaged bone surface is visible in the 
case of postmortem damage. We excluded stress or accident-related 
fractures. In case of the cranium, we only include traumas that were 
located above the hat-brim line—the line where an imagined hat brim 
touches the head, thatis, the frontal or parietal bone—which is usually 
the part of the skull hit by violent blows, in contrast to accidents that 
usually affect the area below the hat-brim line (for more details about 
the accident-versus-violence differentiation, see Supplementary Meth- 
ods). Still, we need to acknowledge that it is not possible to distinguish 
accidental and violent trauma perfectly. We included only skeletons 
aged 15 years and above to obtain a more homogeneous sample based 
on the age criteria. In addition, we collected data about postcranial 
marks typically resulting from weapon wounds, such as arrowheads 
(Supplementary Methods). This makes our study compatible with 
the Global History of Health project studies”. However, the share of 
postcranial weapon wounds was tiny and most of our evidence comes 
from cranial trauma. Although cranial and weapon traumas are obvious 
indicators of violent behaviours and allow us to measure the intensity 
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Fig. 4 | The estimates of the number of buried weapons between the 
Chalcolithic and the Middle Bronze Age. For this figure, we divided the 
estimated number of weapons by summed population (estimated on the basis of 
excavation site areas). Weapon estimates were adapted from ref. 92, population: 
see Assessing the representativeness (or selectivity) of the violence dataset. Error 
bars cannot be provided as we do not have access to the individual observations. 


of interpersonal violence in our sample, it must be stated that such 
a strategy typically underestimates violent deaths since a variety of 
wounds, particularly those to the soft tissues (for example, poison, 
septic wounds and ruptured organs), do not leave specific marks on 
bones. Notably, ref. 78 demonstrated that these data agree well with 
other proxies for interpersonal violence, suchas historical records on 
homicide, regicide and elite deaths in battle—the last include violence 
against soft tissues of course, suchas stabbing a knife into the stomach. 
We calculated trauma shares, which correspond to the ratio between 
the number of cranial and postcranial weapon wounds and the total 
number of sufficiently well-preserved skeletons from a given site/ 
period. Our sample consists of 3,539 total individual human remains 
examined with an overall number of traumas attested on 323 indivi- 
duals, a weighted average share of 9.1% (Table 1). 

Turning to potential selectivity biases that may affect our trends, 
we controlled for the distribution of observations between urban and 
non-urban (burial sites, rural villages, pre-urban communities and so 
on) environments. We created dummy variables for special contexts, 
such as sacrifice of human retainers in royal interments. These vari- 
ables allow us to adjust for religious and royal samples in the following 
regression-based estimates. Moreover, we identified clearly defined 
battle sites or cases of military attacks. In case of violence occurrence, 
the division between urban and rural contexts is considered meaningful 
by a substantial strand of historical literature”’”’. Overall, our sample 
is well-balanced with 1,310 cases coming from urban contexts and 
2,229 from non-urban ones, with this latter category grouping rural 
communities, burial sites and pre-urban communities. 


Assessing the representativeness (or selectivity) of the 
violence dataset 
We assess, in this part of the study, to which degree the selectivity of 
burials changed over time. We cannot claim full representativeness, as 
probably no bioarchaeological sample for sucha comprehensive (early) 
period and region can be fully representative. In the following, we dis- 
cuss at which points our sample deviates from full representativeness 
and how substantial these differences are. Fortunately, we can provide 
alarge body of evidence and plausibility assessments that indicate that 
our sample does not fundamentally deviate from this ideal. 

We could imagine a bias because in some societies individuals 
who were more at risk of violence were better represented among the 
burials that we can observe. For example, in societies that emphasized 
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and idealized military values, their soldiers might have been buried 
with more care and perhaps been more likely to be found by modern 
archaeologists. We call this a ‘potential warrior idealization bias’. How 
can we assess whether such a bias existed and was substantial? Several 
possibilities exist. First, one strategy is to discuss whether, during some 
periods and in regions of our study area, military success was exception- 
ally highly considered, even idealized, and whether this resulted inthe 
practice of burying military individuals with greater care. This might 
lead to more individuals with violent trauma among the overall number 
of buried persons recorded in excavation reports. In fact, some qualita- 
tive changes in military presence took place in society, especially when 
metal weapons became more widespread in graves during the Bronze 
Age and graves of warriors became richer and burials were performed 
with greater care. The evidence strongly indicates that, although swords 
and metal weapons were already introduced during the Chalcolithic 
period (when feuding turned into warfare®*°), weapons start to be 
routinely found in the graves during the EBA-MBA. Several leading 
scholars elaborated on the creation of an elite warlike ethos as reflected 
in burials” ®. Bronze weapons gave warriors a special status”. 

Gemez surveyed all the bronze weapons excavated in the Ancient 
Middle East®’”*. Most of the weapons were excavated in warriors’ burials 
and their frequency takes off after 3,000 BCE (Fig. 4). Notably, this 
increase is dramatically larger than the increase in overall population. 
For the pre-Chalcolithic period, we need to adda cautionary note: many 
weapons were tools that were mostly used for other purposes such as 
hunting, farming and so on but that could also be used as weapons. At 
the same time, during most of the Bronze Age we observe a relatively 
modest and declining violence level, so this bias cannot make our 
findings obsolete; Bronze Age violence levels might be even lower if 
we take this high status of warriors into account. In other words, we see 
an inverse trend in our trauma data. Exactly when the warrior burials 
were richer and more elaborate, which Gernezinterprets as reflecting 
asociety with strong warrior idealization, cranial traumas and weapon 
wound evidence decreased, suggesting less interpersonal violence. In 
sum, we find arguments that the warrior idealization bias was probably 
not too substantial, as for the most extreme case it did not hold. We 
could even put forward the argument that the bias would reinforce 
our results. Moreover, we control for battle sites and massacres, so the 
bias might not be large. Of course, this corresponds to our discussion 
about substituting interpersonal violence with violence against foreign 
states. The Bronze Age states succeeded in reducing interpersonal 
violence but they channelled the aggression of the military parts of 
their societies against other states and societies”. Hence, another 
possible interpretation would be the specialization of warrior activi- 
ties, which leads to less violence outside the warrior occupation™”. 

Next, we pursue a second strategy to assess potential bias by com- 
paring different regions (Fig. 2). The cross-sectional differences of 
potential excavation bias can be assessed by comparing similar trends 
for different regions. In our case, we observe graphically a strong 
inter-regional co-movement of violence rates and it seems less likely 
that all these societies had similar levels of warrior idealization bias. 
While the changing selectivity could be imagined having taken place 
in these other regions as well, such synchronicity is less likely under 
different institutional and social circumstances. Hence, similar trends 
increase the plausibility of individual trend estimates. 

Thirdly, the correlation with possible determinants of interper- 
sonal violence makes the estimates even more plausible. The literature 
onthe Ancient Middle East clarifies that the EBA and MBA were periods 
in which the states had developed skills to convince their people that 
less interpersonal violence should be used, as, for example, increas- 
ingly codified legal systems helped to solve conflicts without violence” 
(Supplementary Appendix A). As a result of these new capabilities of 
states, inter-regional trade flourished. This correspondence of poten- 
tial causes and consequences makes it substantially less likely that our 
observed trends stem from social selectivity of excavated individuals. 


All three strategies of bias assessment add additional plausibility to 
the newestimates. 


Empirical approach to deteremine variables related to 
violence trends 

In the following, we formally assess in a multiple regression analysis 
whether the changes over time that we observed were statistically sig- 
nificant (see notes to Table 2). Moreover, we include control variables 
to analyse whether the changes towards less interpersonal violence 
during the EBA-MBA might have been correlated with third variables 
such as climatic shocks, population pressure or military innovations, 
instead of the correlation with the period variables EBA-MBA. Follow- 
ing the extant research on violence, it is possible to single out many 
factors as potential covariates, such as state capacity, demography, 
income, inequality, trade, human capital, military technologies, cli- 
mate, diet, geography and so on””*”’. We focus here on some key fac- 
tors that can be reliably measured in the archaeological record and we 
evaluate whether our estimates are affected by them or not. We control 
for settlement size”, the most important military innovations (iron 
technologies and mounted warfare’), climate shocks (droughts, that 
is unusually cold/dry spells”) and population fluctuations. Our main 
control variables are defined as follows: ‘urban versus rural’ character 
is Lifa threshold of 5 ha area size is reached and O otherwise. ‘Climatic 
shocks’ are defined as the share of the period that is characterized 
by a rapid climate change event consisting of cool and dry spells”, 
reducing temperature and precipitations and generating droughts 
that negatively affected agricultural production”. 

How do we measure weather shocks? Several studies have focused 
on palaeoclimatic evidence, especially in the debate about the trans- 
formation towards agriculture and urbanization, in the debate about 
the collapse events 4.2 and 3.2 ka before present and on the impact 
of climate change and state evolution. Palaeoclimatic evidence 
routinely exploited by this literature consists of natural archives of 
high-resolution proxies (speleothems) or global simulated trends 
scaled at various spatial resolutions. According to growing evidence 
about the positive correlation of lethal violence and drought-induced 
resource shortages in the pre-industrial world"”*?> in our study, we 
focus oncontrolling for the impact of droughts (that is, unusually cold 
and dry climatic spans) which were negatively related with production 
by lengthening the dormant phase of plant growth (winter crops in 
the study sample), creating mismatches in the agricultural cycle”. 
Although local conditions may vary slightly, extended droughts (such 
as millennial-scale events) tend to be correlated between regions. 
Accordingly, we wanted to be able to check whether, in the aftermath 
of major drought spells, violence spiked in our sample. 

‘Iron technology’ is defined as 1 if the excavated site yielded iron 
materials and especially iron weapons. This military technology was 
introduced first in the Caucasus region and then imported into Meso- 
potamia® and other regions of the Middle East. It might have been 
related to social change in our regions. Similarly, we add a control vari- 
able for the emergence of horse-mounted warfare. Horse-based mili- 
tary technology arrived from Central Asia and the Caucasus region”. 
This variable is 1 if horses are shown in pictorial finds, mentioned in 
cuneiform tablets or bioarchaeologically attested”. Finally, ‘popula- 
tion size’ is included in the equation. We use the estimates of population 
size based on the area of the excavation locations during a specific 
time period, to control population pressure effects. We do not make 
any claims about the direction of causality or the exogenous nature of 
these variables. We only want to assess whether these variables might 
have made the change from the Chalcolithic to the EBA-MBA statisti- 
cally insignificant. 


Main results of regression analysis 
In Table 2 we provide the regression results, for which we have added 
control variables and assessed statistical significance. The reference 
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category (the comparison period) for the time periods is the Chalco- 
lithic, a period with tremendous violence as seen in Fig. 2. Here, we 
briefly summarize the results again to compare them with models in 
which we are adding control variables. In comparison with the Chalco- 
lithic period, we observe lower violence during the Neolithic period 
(Table 2, column 3: b = -1.332, z(70) = -4.33, P= 0.000, Cohen’s d effect 
size 1.193, 95% CI -1.934, -0.729). Again, relative to the reference cate- 
gories of the Chalcolithic, we observe a decline in violence during the 
EBA and even more pronounced in the MBA (Early Bronze b = -0.800, 
2(70) = -3.22, P= 0.001, Cohen’s deffect size 0.741, 95% CI -1.287, -0.312. 
Coefficient b = -1.815, test statistic z(70) = -2.55, P= 0.011, Cohen’s d 
effect size 2.852, 95% CI -3.210, —0.420). In contrast, while for the IA 
the coefficient is also negative (hence less interpersonal violence than 
inthe Chalcolithic), the coefficient is smaller (b = -1.114, z(70) = -4.38, 
P=0.000, Cohen’s d effect size 0.997, 95% C1 -1.612, -0.615). The coef- 
ficients are reported in logs. If we want to express the effect in percent- 
age, we need to take the exponential. For example, the coefficient of 
MBA (relative to Chalcolithic) is exponential(1.815) = 6.141, witha nega- 
tive sign. Hence the MBA experienced a decline in violence by 6.141%, 
whereas the Iron Age only saw lower level of violence compared to the 
Chalcolithic. To make sure, we tested the Holm-Bonferroni adjustment 
in Supplementary Appendix C and arrived at the same results as the 
results reported here, without adjustment. 

We tested a model without controls (Table 2, column 1) and then 
added the controls one by one. The results of the EBA and MBA period 
coefficients turn out to be extraordinarily robust. In model 2 of Table 2, 
we added regional differences (regional fixed effects) to control for 
differences between the four regions that could not be covered by other 
proxy variables. In column 3 of Table 2, we add variables to control for 
a battle character of the site, a royal interment or an urban character. 
The royal and battle variables are not tested separately in separate 
columns because we do not include them with a contextual model in 
mind but simply as controls. Royal and battle variables are associated 
with higher violence (b = 1.157, z(70) = 2.600, P= 0.009, Cohen’s deffect 
size 1.134, 95% CI 0.283, 2.031) and royal context (b = 1.861, z(70) = 3.970, 
P=0.000, Cohen’s d effect size 1.824, 95% CI 0.941, 2.781). In column 
4 of Table 2, we add climate shocks, in columns 5 and 6 the iron and 
horse weapon technologies and in column 7 population size. Adding 
the control variables one by one, we observe that the coefficients for 
EBA-MBA are almost unchanged. Urban sites had no substantially 
higher interpersonal violence during this period. In column 4 of Table 2, 
we add the variable for climate shocks. The coefficient of this variable 
is statistically insignificant. In Table 2, columns 5 and 6, we add iron 
and horse-riding technology. Again, the EBA-MBA coefficient was not 
substantially affected. 

To assess the robustness and sensitivity of our results, we perform 
a series of tests in Supplementary robustness and sensitivity tests: 
(1) we begin with an assessment of time-region units with small num- 
bers of observations; (2) we also take a different modelling strategy for 
our multiple regression analysis—we do not cumulate the right-hand 
side variables but include them separately; (3) we performa sensitivity 
analysis by including an interaction term analysis of urban place and 
period. This allows us to assess the robustness of our results consider- 
ing urban and rural differences (or not). We finish with (4) an anthro- 
pological sensitivity analysis, namely, including or not including the 
so-called parry fractures (Supplementary Methods). We conclude from 
all these robustness and sensitivity analyses that our results are quite 
stable and consistent. 

In sum, the observed coefficients correspond with the propos- 
als that we formulated on the basis of the existing literature and the 
additional control variables did not change the results. During the 
Neolithic period, when population density was very low and therefore 
potentially conflicting groups could easily avoid each other, interper- 
sonal violence was low. Moreover, we formulated the proposal that 
during the EBA-MBA the rule of law and the synergetic effect of law with 


trade expanded, implying also relatively low interpersonal violence. 
Again, this was confirmed, even after controlling several other vari- 
ables. Although we cannot generalize these results for other regions 
and periods, we can confirm that the proposal corresponds with the 
observed development of violence in the pre-Classical Middle East. 


Reporting summary 
Further information on research design is available in the Nature Port- 
folio Reporting Summary linked to this article. 


Data availability 

The data that support the findings of this study are available in the 
publicly accessible Zenodo repository at https://doi.org/10.5281/ 
zenodo.8010025 and inthe Supplementary Data. 


Code availability 

Custom code that supports the findings of this study is available in 
the publicly accessible Zenodo repository at https://doi.org/10.5281/ 
zenodo.8010025. The Stata code is also available in the Supplementary 
Code. 
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